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Abstract 

The radial part of the Klein-Gordon equation for the generalized Woods-Saxon potential is 
solved by using the Nikiforov-Uvarov method in the case of spatially dependent mass within the new 
approximation scheme to the centrifugal potential term. The energy eigenvalues and corresponding 
normalized eigenfunctions are computed. The solutions in the case of constant mass are also studied 
to check out the consistency of our new approximation scheme. 
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I. INTRODUCTION 



The investigation of the quantum mechanical systems in the case of position dependent 
mass (PDM) following of works by v. Roos, and Levy-Leblond [1, 2] have recently been 
received great attentions. This is so because such solutions are aviable in wide range of 
different areas, for example, in the study of impurities in crystals [3-5], or of electronic 
properties of quantum wells, and quantum dots [6], and in semiconductor hctcrostructurcs 
[7]. Yahiaowi, and Bentaiba [8] have studied the weak-pseudo-Hermiticity in the case of 
PDM, Ganguly, and Nieto have extended the second-order supersymmetric approach to 
the systems with coordinate dependence mass [9]. In Ref. [10], some new shape- invariant, 
exactly solvable potentials are generated by using a specific ansatz in the point of PDM-case. 
Ju et al. have been studied the dynamics of a quasi-free particle in an effective potential 
arising from the dependence of the mass on coordinates, and analyzed the eigenfunctions 
and probability densities for s- waves [11]. 

Another interesting area received a lot of attentions is that solving the Schrodinger (SE), 
and Dirac equations in the case of PDM. To solve the above equations has been used different 
methods, and approaches for different potentials, such as deformed algebras in Coulomb po- 
tential [12], in the content of supersymmetric quantum mechanics [13, 14, 15, 16], quadratic 
algebra approach [17], numerical analysis of a square potential by using appropriate match- 
ing conditions [18], point canonical transformation applying on harmonic oscillator. Coulomb 
and Morse class of potentials [19], finding the non-relativistic Green's functions with PDM 
for harmonic oscillator [20], Coulomb potential in Dirac equation [21], Morse potential in 
PDM background [22] , a series solution of the SE for Cornell potential [23] , finding the bound 
states of Rosen- Morse and Scarf potentials via general point canonical transformation [24] . 

In the present work, wc give the approximate solutions, and corresponding wave functions 
of the radial Klein-Gordon (KG) equation for the Woods-Saxon (WS) potential in the case 
of PDM. We investigate the energy spectrum, and the corresponding eigenfunctions of the 
generalized WS potential by using a new approximation to the centrifugal potential. In order 
to find the spectrum we use the NU-method in the case of exponentially mass distribution 
varying with coordinate. The NU-method is a powerful tool to solve of the second order 
linear differential equations with special orthogonal functions. In this method, the differen- 
tial equation is turned into a hypergeometric type equation by using a transformation on 
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coordinate [25]. 

The organization of this work is as follows. In Section II, we solve the radial part of 
the KG-equation for generalized WS potential by using the NU-method within the frame- 
work of an approximation to the centrifugal term. We compute the energy eigenvalues and 
corresponding eigenfunctions, and also give the results for the case of the constant mass to 
control the consistency of our new approximation. We write our conclusions in Section III. 



II. NIKIFOROV-UVAROV METHOD AND CALCULATIONS 

In spherical coordinates, the radial part of the Klein-Gordon equation can be written as 



where i is the angular-momentum quantum number, E is the energy of the particle, m is 
the rest mass, and c is the velocity of the light. 

The generalized WS potential can be written of the form [31] 



V{r\ = ^- , (0 < r < oo) . (2) 

where Vq is the potential depth, /3 is a short notation, i.e. /3 = 1/a, a is diffuseness of the 
nuclear surface, r is the ccntcr-of-mass distance between the projectile and target nucleus, 
and To is the width of the potential, which is proportional with target mass number A. q is 
the deformation parameter, and arbitrarily taken to be a real constant. The WS potential is 
widely used in the coupled-channels calculations in heavy-ion physics. This model explains 
the single-particle motion during a heavy-ion collisions [26-29] . 
Let us write the potential as 



V{x) = - V ' (3) 

where x = (r — ro). Eq. (1) can not be solved exactly because of the centrifugal potential 
term iox I . The nuclear distance r can not fluctuate very far from the equilibrium for 
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rather high vibrational levels [32], which gives small x- values. So the centrifugal potential 
term can be expand about x — ets the following 



£(£+1) ^^^^ X. 



+ 3 



X 



To 



+ ...), 



(4) 



where the parameter D in the above equation is given as D 



2mri 



Instead, we suggest to replace Vi{r) by the following potential form [33] 



(5) 



where the parameters Dq,Di , and D2 are arbitrary constants. 

Expanding the potential V-[{x) around x = under the same condition, and than com- 
bining equal powers with Eq. (4), one can find the arbitrary constants Di{i = 0, 1, 2) in the 
new form of the potential as 
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(6) 
(7) 
(8) 



where it can be seen that the new parameters Do,Di , and D2 are real, dimensionless param- 
eters, and dependent to the numerical values of the quantum system under consideration. 
On the other hand, we prefer to use the following position dependent mass function 



m{x) = mo 



1 _ (1 + ^e'^^ 
mo ^ 



(mo > mi) , 



(9) 



where mo and nii are two arbitrary positive parameters. The mass function is finite at 
infinity, and enables us to solve analytically the KG-equation given by Eq. (1), and to check 
out the limit of the case of the constant mass. 

Substituting Eq. (5), and (9) into Eq. (1), we get 
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h^i{i + 1) 



Dn 



1 



(E'H + moTTiic ) - 



h^i{i + 1) 



1) 



1 + ge/^=^ 
(l + Je/^.)2}'^(^) = 



By using the transformation z — 2(1 + qe ) , we have 



(f^{z) 2(1 - z) d(f){z 



z{2 - z) dz [z{2 - z)\' 



a\z^ — a\z — at 



where 



u;le{e + 1)D2 + uJi{mY - V?) , 
2\uj\t[l + \)D^ - 2cj2(momiC^ + EVq)\ , 
^{u)\t{t + 1)^0 + - E")] . 



and Jl = 1/I3^rl , and = . 

To apply the NU-method, we rewrite Eq. (11) in the following form 



(10) 



(11) 



(12) 



0"(^)+44</''(^)+4^0W-O, 



(13) 



a{z) cr^(-2) 

where (t{z) and alz) are polynomials with second-degree, at most, and f{z) is a polynomial 
with first-degree. By using the following transformation for the total wave function 



<P{z) = az)i^{z) 



(14) 



we get a hypergeometric type equation 



a{z)ilj"{z) + r{z)ilj'{z) + Xilj{z) = 0, 
where ^{z) satisfies the equation 
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(15) 



i\z)/i{z)^7:{z)/a{z). (16) 

and the other part, 'ip{z), is the hypergeometric type function whose polynomial solutions 
are given by 



where 6„ is a normalization constant, and the weight function p{z) must satisfy the condition 



j-J,a{z)p{z)]=T{z)p{z). (18) 
The function t:{z) and the parameter A required for this method are defined as follows 



, , a'{z)—T{z) , ,a'(z) — t(z) , , . s ..^s 



A = A; + Ti\z) (20) 

The constant k is determined by imposing a condition such that the discriminant under the 
square root should be zero. Thus one gets a new eigenvalue equation 



A = K^-nr'- (^^0,1,2,...) (21) 

where 



t{z) = f{z) + 2t:{z) . (22) 

and the derivative of t{z) must be negative. 
Comparing Eq. (11) with Eq. (13), we have 

f{z) = 2(1 - z) , a{z) = z{2 - z) , a{z) = -a\z^ - a\z - a\ (23) 

Substituting this into Eq. (19), we get 
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7r{z) = ±^ {aj - k)z'^ + {al + 2k)z + a 



(24) 



The constant k can be determined by the condition that the discriminant of the expression 
under the square root has to be zero 



(a| + 2kY - Aai{ai - k) ^ 

The roots of k are ki^2 = ~ ^ (^2 ~ ^ ^ dsA, where A 
these values into Eq.(19), we get for 7r{z) for ki 



(25) 



al + 2a2 + 4ai. Substituting 
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(26) 



and for /c2 



7r(^) = T 



/ ^ , "3 \ 
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Now we find the polynomial t{z) from ti{z) for the second choice as 



(27) 



T(z) = 2 + 2a3-2(| + ^ + l)z. 
so its derivative — 2^^ + ^ + l) negative. We have from Eq. (20) 



(28) 



^ = - ^ {al + al + Aa3 + A + as) , 



(29) 



and Eq. (21) gives us 



A» = 2n(4 + |+l) + 



n — n. 



Substituting the values of the parameters given by Eq. (12), and setting A 
find the energy eigenvalues for any ^-states 



(30) 



A„, one can 
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n,i 



Vq[4u;^V^ + N^- Aul£{£ + l)(Di + D2) + AmiC^u^{2mo - mi)] 



2{m + Au}iVi) 



U>2 



ul{2Do + Di + D2W + 1) + 2mgcM 



'ul{Di + D2W+l) 



n 2 



+ /^ojlVi 



(31) 



where the energy eigenvalues with (+) sign correspond to particle, and the one with (-) sign 
correspond to antiparticle. Two parameters in the above expression are 



mi 



and 



8mic4u;i(2mo - mi)[4a;i2(L>i + D)^{^ + 1) - 2ujlmic*{2mQ - mi) - {N^ + AulVi) 



(32) 



N = -{2n + l) + Jl + Aa\. 



(33) 



We see that the energy levels for particles and antiparticles are symmetric, and the ground 
state energy is different from zero. We summarize some numerical results in Table I to see 
the effect of the spatially dependent mass parameter mi on the energy eigenvalue of bound 
states. It is observed that the energy levels are strongly dependent on the parameter, and 
the increase of the energy eigenvalues in the existence of mi is very significant. It has to 
be stress that the higher numerical values of the parameter mi give positive values for the 
bound states. 

The energy spectra in the case of constant mass is obtained by setting mi = in Eq. 
(32) which gives us mi = , and we get 



e: 



mi=0 



Vq[AuIV^ + iV2 - AujII{1 + l)(L'i + D^)] 



2{N^ + AujlVi) 



± 



0J2 



\ 



u;l{2Do + Di + D2W + 1) + '^rnlc^ul 



' u;f{Di + D2)e{e + 1) 



n 2 



where 



N' = -{2n + 1) + \ll + 4af , af = a?(mi ^ 0). 



1 

16 ' 
(34) 



(35) 
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It is seen that the result for the case of constant mass is the same with those obtained in 
Ref. (29). 

In order to find the eigenf unctions, we first compute the weight function from Eq. (18) 



p{z)^z^^{2-z)^, (36) 

and the wave function becomes 

where bn is a normahzation constant. The polynomial solutions can be written in terms of 
the Jacobi polynomials [34, 35] 

V^„a^) = ^'ni"i"^'^^(l-^), A>-1, a3>-l. (38) 
On the other hand, the other part of the wave function is obtained from Eq. (16) as 

^{z) = ^"3/2 (2 - z)^/^ . (39) 
Thus, the total eigenfunctions take 

<Pr^e (Z) =b'^{2- ;,)^/2^»3/2 p(as, A) _ ^) _ (4q) 

where b'^ is the new normalization constant. It is obtained from 

f [' \Mz)f {'-^)dz ^ 1 . (41) 
To evaluate the integral, we use the following representation of the Jacobi polynomials 



n\r{n + a + <; + !) 

" (n\ r{n + a + c; + m + l)r{n + a + <; + m + l) z-l 
i^,[J r{m + a+l) r{m + a + l) ^ 2 > ' ^ > 
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where j ^ r\{n-ry. ^ r(r-+i)r(t-r+i) ■ Hence, from Eq. (41), and with the help of Eq. 
(42), we get 



[g(n, m) x ^(r, s)][^) \h'S [ ^-+^+'^3-1 (2 _ ^)^+i dz^l, (43) 



where g{n, m), and g{r, s) are two arbitrary functions of the parameters A, and 03, and given 

by 



2-™r(A + n + l) 

g[n,m) 



n!r(A + a3 + n + l) 

/ n \ r(n + a + ^ + m + 1) r(A + 03 + n + m + 1) 
""ioUy r(m + a + l) r(a3 + m + l) ^ ^' 

(44) 

and 

5f(r, s) = g{n, m) {n ^ r;m ^ s) . (45) 

The integral in Eq. (43) can be evaluated by using the following integral representation 
of hypergeometric type function 2-^1(0, b; c; z) [36] 

2F,{a,b;c;z)= {1 - tr'-' {1 - tz)''^ dt , (46) 

by setting the variable z ^ |, and taking c— l + b,z — 1, one gets 

|*'-(2-*r*=^^^.F.(a,M + ^;i), (47) 

From last equation 

^m+.+a.-l /o „nA+1 j„ r(m + S + a3)r(l) 



2«r(m + s + a3 + 1) 



1 



X 2-Fi(-^ - 1, m + s + 03; m + s + 03 + 1; - ) , (48) 
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where we set b — m + s + , and a — —A — 1. 

By using the following identities of hypergeometric type functions [36] 



r ife + i r(6-a + i) 

2Fi(a,6;c;-l) = J' , (a - 6 + c = 1 , b > 0) (49) 

2F1 (a, 6; c; ^) = 2" ai^i (a, c - 6; c; -1) , (50) 
the function of 2-^1 (a, i*; c; 2;) in Eq. (48) can be evaluated as 



2F1 {-A - 1, 6; 1 + 6; 1) = 2-^-' 7f ^ !l ' (m + 5 + a3-^ = 2), (51) 
Finally, we get the normalization constant as 



^ ^ r(»^ + .^ + »:. + i)r(f +A) 

' /30Fr(m + s + a3)r(3 + A)[y'(n,"i) x^(^,s)] ■ 

where 



^'(n, m) 



T{A 



n 
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n!r(A 



X 




+ 0,3 + n + 1) 

n \ T{n + (T + ^ + m + l) r(A + 03 + n + m + 1) 
r / 



r(m + a+ 1) 



r(a3 + m + 1) 



■1) 



m+l 



(53) 



III. CONCLUSION 



We have solved the radial part of the KG-equation for the modified Woods-Saxon po- 
tential in the case of position dependent mass by using a new approximation scheme to the 
centrifugal potential term for any i values. It is observed that the results obtained by using 
the new scheme for the case of the constant mass arc consistent with the ones obtained in 
Ref. [29]. It is seen that there is a linear relation between the energy eigenvalues and the 
contributions coming from the dependence of the mass on spherical coordinate. The energy 
spectra and the corresponding wave functions are obtained by applying the NU-method. 
The eigenfunctions can be expressed in terms of Jacobi polynomials in the scheme of the 
new approximation to the centrifugal barrier in the case of position dependent mass. 
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TABLE I: The dependence of the bound states for a system 'proton+nucleon with 


average mass 


number A = 56' on the parameter nii in AleV for q = 


1 by using the numerical values ra^ = 


1.007825 amu, Vq = 47.78 MeV, ro = 4.91623 fm [37]. 
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